Abstract 



The main topics of this paper is to shown a Fast Mie Algorithm 
FMA as the best way to use the Mie scattering theory for cross section 
calculation. This fast algorithm used recursion for summing a long 
timed sum of cylindrical functions. 
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1 Mie Theory 



Gustaw Mie in 1908 had considered a problem of the light scattering by the 
macroscopic dielectric spheres. Mie find an exact solution of the problem of 
diffraction the plane wave on the homogenous sphere with constant density. 
Similar solution was published by Debay. 

The Mie solution satisfies also the scattering on many spheres with the 
same structure and diameter - like scattering on the drops. The Mie solutions 
in widely used in various situations up to nowadays. 

Many scientists (0, [@], @, [0 and others) show a lot of simply models, 
more simply than Mie 's to application. Unfortunately, all of this approxima- 
tions, works only in strictly defined area (sometimes with great error). For 
precisely calculation done in all the spectrum of measure we must used the 
main theory, now is only one such theory - Mie Theory. 

On the long distances, the process of light scattering is described by 
amplitude A of electric field E: 

as exp(ikr) 

E = A— -, kr » I. (1) 

ikr 

The dependencies of angle give us angle distribution like: 

da/dVt = \A/ik\ 2 . (2) 

After integration over all angles it takes us total cross section. Amplitude 
can be write also in such form: 

A = nx(exn)ff!+nx(exk)ff 2 (3) 

This form is more useful for calculation of physically dependencies. For 
example, amplitude of radiation, with parallel polarisation, is equal to: 

A\\ = (H 1 cos 6 + H 2 ) cos (p, (4) 

and for perpendicular 

A x = - (Hx + H 2 cos 9) sin ip . (5) 
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For non-polarised wave the intensity of reflected light is equal: 



,2 



h(0) = ~\Hicos6 + H 2l 
I±{6) = ^|#i + # 2 cos#| 2 
These components give us total intensity 1(9) : 



1(0) = 1 + c ^ e (1^(2 + \H 2 \ 2 )+2cos6Re(H 1 HZ) (6) 



2 

Polarisation P(#) is given as: 



|#i| 2 - \H 2 \ 2 



If i/ 2 = 0, then we get the known Rayleigh's results for amplitude 

A R = n x (e x n)ifi. (8) 

Intensity is expressed by: 

/fi(S) = i±|^|^, (9 ) 

and polarisation as: 



2 The Mie sum for amplitude H 

We consider electrical field satisfied the Maxwell equations ||, |J, [0, f|. If 
£ = cos-#, then: 



Hi = (M — tE)f, H 2 = {E-tM)i. (11) 
Using an operator notation |7J we can write: 

Di p) = 6 p h x (a) 3i (fi) (12) 

M = fidp - ad a 

Oe = ocd/s — fid a 

a = kR 

fi = NkR. 

Both E and M have a well-known structure: 



00 9/_i_i 

"=Ei(fV*w*. (14) 

where P; (i) are the Legendre polynomials. Now, E and M can be written as 

E l = e«?sin6?=*^, (15) 



;m . M ReZ) 7 A/ 

Gin _ — 



M^e^r sin 5 f = _ j*-, (16) 

where ji(z), hi(z) are cylindrical functions of Bessel and Hankel (first kind 
functions). We can write now the electric part 

Jl (a)f l (fi)-Nf l (a)j l (fi) 
1 hiaW) - Nh\(aWy {i> 



and magnetic one 



N3i{a)jm-3[{a)m) (m 
1 NhMjUfi) - h'iaWY 1 ] 
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3 Approximations 



For small a, the scattering cross-section was written by Lord Rayleigh |L 

pi, m 



Q 



8a 4 (N 2 - l) x 
3 (N 2 + 2f 



(19) 



For water, this approximation work for a < 0.5 with error smaller than 1%. 
In range of a » 1, good approximation take us the formula of van de Hulst 



Qh 



sin 5 1 — cos 5 
1 - 2—^- + 2 



S ■ 5 
Scattering cross - section is defined as: 



S = 2a(N-l) . 



(20) 



a sc = I d{l\A/k\ 2 . (21) 



After normalisation by ira 2 is equal: 



9 , . 

g = ^^(2/ + l)(|^| 2 + |M i | 2 ). (22) 
a i=i 

So, we obtained total cross section. 



4 FMA - Fast Mie Algorithm 

Our program calculated Eq. (P2"D . This program is written in Turbo Pascal 
language. The structure of algorithm is very simple so the translation to 
another language can take no more like some minutes. 

In program we used recursion, so the time of calculation is minimized 
to some seconds for thousand of dates. Recursion used main properties of 
calculated functions and is most efficient for great alpha. 

Program generates the data file with location C:\MIE.DAT. 
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PROGRAM FASTMIE; {simplified source code} 

Uses Crt; 

Var 

Wy:Text; 

wave_k, wave: Integer; 

q,n,w,alpha,alpha_k,alpha_p,alpha_step : Extended; 

crossection, crossection_el, crossection_ml :Extended; 

a,b,c,d,e,f, pp, m :Extended; 

jO, jl, j2, jbO, jbl, jb2, yO, yl, y2 Extended; 

Begin 

ClrScr; 

Writeln ('Program FAST MIE'); 
Writeln; 

Writeln ('Give N (water N=3/4)'); Readln(n); 
Writeln ('From alpha '); Readln(alpha_p); 
Writeln (To alpha '); Readln(alpha_k); 
Writeln ('Step of alpha '); Readln(alpha_step); 
Assign (Wy,'C:\MIE.DAT') ; 
Rewrite (Wy); 
alpha:=alpha_p; 
Repeat 

alpha :=alpha+alpha_step; 
crossection :=0; 

{ calculation for 1 st wave } 
wave_k:=round (alpha* 1.12+10); 

{initial values of cylindrical functions } 

{ j=Re_h y=lm_h } 

jO:=sin(alpha); yO:=-cos(alpha); 

jl:=jO/alpha+yO;yl:=yO/alpha-jO; 

jbO:=sin(alpha*n);jbl:=jbO/(alpha*n)-cos(alpha*n); 

pp:=l; 

a:=yl; b:=y0; 
c:=jl; d:=j0; 
e:=jb0; f:=jbl; 
w:= (n-(l/n))/alpha; 
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{electric part } 

m:=((a*e)-(n*b*f)+(w*a*f))/((c*e)- (n*d*f)+(w*c*f)); 
q:=(2/(alpha*alpha))* (pp+2); 
crossection_el:=q*(l/(l+(m*m))); 
crossection:=crossection+crossection_el; 
{magnetic part } 

m:=((n*a*e)-(b*f))/((n*c*e)-(d*f)); 
crossection_ml:=q*(l/(l+(m*m))); 
crossection:=crossection+crossection_ml; 

{recursion for all partial waves } 
For wave: =2 to wave_k Do 
Begin 

pp:=(2*wave-l); 

j2:= pp * jl / alpha - jO; 

y2:= pp * yl / alpha - yO; 

JO: jl: jl: j2; 
yO:=yl; yl:=y2; 
a:=yl;b:=yO; 
c: =jl;d:=jO; 

jb2:=pp*jbl/(alpha*n)-jb0; 
jbO:=jbl; jbl:=jb2; 
e: =jbO; f:=jbl; 
w:= wave*(n-(l/n))/alpha; 
{electric part } 

m:=((a*e)-(n*b*f)+(w*a*f))/((c*e)-(n*d*f)+(w*c*f)); 
q:=(2/(alpha*alpha))* (pp+2); 
crossection_el:=q*(l/(l+(m*m))); 
crossection:=crossection+crossection_el; 
{ magnetic part } 

m:=((n*a*e)-(b*f))/((n*c*e)-(d*f)); 
crossection_ml:=q*(l/(l+(m*m))); 
crossection:=crossection+crossection_ml; 
End; 

Writeln(Wy, alpha:5:5, ' ' , crossection:5:5); 
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Until (alpha>=alpha_k); 

Close(Wy); 

End. 



5 Results 

Now, we can show some results as plots for N equal to 4/3, 1.1. On all these 
plots we can see very interesting structure with resonances Fig 3,4 and 5. 

Fig. 1 The cross - section for N=l.l for alpha from to 100. 
Fig. 2 The cross -section for N=4/3 for alpha from to 100. 
Fig. 3 The cross - section for N=4/3 for alpha from 60 to 65. 
Fig. 4 The cross-section for N=4/3 for alpha from 63 to 64. 
Fig. 5 The cross-section for N=4/3 for alpha from 63.30 to 63.34. We 
can see the Lorenz curve. 
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